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The occurrence of single- or multisite localized vibrational modes, also called Discrete Breathers, 
in 2D hexagonal Dusty Plasma lattices is investigated. The system is described by a Klein-Gordon 
hexagonal lattice characterized by a negative coupling parameter e in account of its inverse dispersive 
behavior. A theoretical analysis is performed in order to establish the possibility of existence of 
single- as well as three-site discrete breathers in such systems. The study is complemented by a 
numerical investigation based on experimentally provided potential forms. This investigation shows 
that a dusty plasma lattice can support single site discrete breathers while three-site in phase 
breathers could exist if specific conditions, about the inter-grain interaction strength, would hold. 
On the other hand, out of phase and vortex three site breathers cannot be supported since they are 
highly unstable. 
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I. INTRODUCTION 



Significant attention has been paid to highly localized time-periodic vibrating motion in periodic lattices since 
roughly two decades ago. These modes [referring to discrete breathers (DBs) and earlier to intrinsic localized modes 
jy-j , (ILMs)], owe their existence to lattice discreteness, in relation to the intrinsic nonlinearity of the medium involved. 

■ Following some initial phenomenological works in the late 80's (see e.g., in Refs. Ql 0, S 0, HI), a rigorous proof 
^\ \ for discrete breather existence was furnished independently by MacKay and Aubry ](| (later extended in Q ) , using 
the notion of the anticontinuous limit, and by Flach [8], who used a homoclinic orbit approach. Although single-site 
breathers possess the lion's share in this study, multi-site excitations (referred to as multibreathers) also exist and they 
are mentioned even since [6| but they were studied more rigorously in later works. Some of the results in existence 
and stability of multibreathers in one dimensional Klein-Gordon chains can be found in Refs. b m el m (also via 
the homoclinic orbit app roach [HI [HI)- For the two-dimensional hexagonal case that interests us one can refer to 
[lH or, for instance, jl6l. Il7l. Il8j for the DNLS case (also, see the references therein). Various studies have therefore 
been dedicated to different aspects involved in th e sp ontaneous formation, mobility and interaction of DBs, both 
theoretically and experimentally; see e.g., in Refs. jl9l. I20I |2ll. I22I I23I |24| for a review. 

On a separate physical playground, large ensembles of char ged particles (plasmas) contaminated by massive heavily 
charged dust defects {dusty plasmas, DP, or complex plasmas) |25ll26ll27ll28| host a wealth of novel linear and nonlinear 
collective effects, which are readily observed in laboratory and space observations. Of particular importance is the 
occurrence of strongly coupled DP crystalline configurations, typically formed and sustained in the sheath region 
above a horizontal negatively biased electrode in gas discharge (plasma) experiments [2H [28l |. These DP lattices 
are encountered in various configurations, including most commonly a superposition of levitated hexagonal two- 
dimensional (2D) layers [28|. One-dimensional (ID) configurations have also been created in the laboratory by making 
use of appropriate confining potentials [2^ | . Ultra- low-frequency modes (eigenfrequencies of only a few dozens of Hz) 
have been predicted, and later experimentally established, in the l ong itudinal (in-plane, acoustic mode), horizontal 
transverse (in-plane) and vertical transverse (off-plane) directions |25l l26l I27I [28| . The linear analysis of the latter 
(transverse, off-plane) degree of freedom has revealed the existence of a backward-wave character, as was theoretically 
predicted [3(|[3l[ and later experimentally confirmed [2§|. Beyond ID "toy-lattice" phenomenology, studies of the same 
off-plane "bending mode" have been carried out recently 2D DP crystals, establishing the inverse dispersive character 
of transverse dust lattice waves |32| . [33|; those findings have successfully interpreted earlier numerical (molecular 
dynamics) and experimental [35| results. It may be pointed out, for rigor, that the 2D picture bears substantially 
modified dispersion and group velocity characteristics, compared to the ID case, as discussed extensively in the latter 
references. The nonlinear modulation of continuum (as opposed to discrete, considered here) dust-lattice excitations 
in hexagonal crystals has been considered for the transverse mode in Ref. [331 ], and for horizontal (in-plane) motion 
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in Ref [3J. 

Here, we aim at focusing on the nonlinear aspects of the transverse (vertical) dust-lattice (TDL) mode, assuming 
the other two degrees of freedom to remain "frozen" throughout this work. The necessary ingredients for discrete 
excitations to occur, namely nonlinearity and lattice discreteness, are present in the transverse dust lattice mode 
par excellence. The sheath environment provides an on-site ("substrate") effective potential (viz. the force is F — 
—d<S>/dz) which may (for low density/pressure) be strongly anharmonic near equilibrium, of the generic form: 

$(z) » $(z ) + ^Mu%(6z n ) 2 + ia (8z n f + {8z n f + 0[{8z n f] . (1) 

The linear TDL eigenfrequency loq [in JT])], is typically as low as m 20 Hz (3(J HH, H3] • The anharmonicity cofficients a 
and (3 may be determined experimentally [38l l39l . |40| • Sources of nonlinearity may include the electrostatic coupling 
anharmonicity [37], EH, l43l. [4a and the coupling among different modes (either due to the interaction law [44| or as 
a purely geometric effect [45j). The anharmonicity (non-parabolic form) and the asymmetry (as in fact manifested 
by a strong cubic term) of the sheath potential form in (1) above is, as a matter of fact, associated with dust-grain 
charging and with the sheath collisionality, as theoretically shown from first principles in Refs. (46[ (note Fig. 8 
therein) and |47| (note Fig. 6 therein). Lattice discreteness, on the other hand, is manifested via the weak interaction 
potential energy among neighboring grains, as compared to the energy "stored" in an isolated (single-site) vibration. 
The former is measured by the characteristic coupling frequency 

V2/1 . *\ V2 



(for Debye interactions) where k = a/ Ad is the lattice parameter (roughly of the order of unity), Q is the electric 
charge of the dust grain and A^> is the Debye length, measuring the strength of Debye charge screening [3(| HH, H3| • 
The latter (single-grain vibration energy) is measured by the linear transverse dust-lattice vibration eigenfrequency 
loq in ([!]). Lattice discreteness is measured via the ratio 

5 = 4o/ w o. ( 3 ) 

which acquires very small values, as suggested by experiments with 2D [S(| and earlier ID [HI, H§] crystals. Dust 
lattices are therefore highly discrete systems, whose dynamics is located near the anticontinuous (uncoupled sites) 
limit (since 8 <C 1). All of these aspects can be inferred from earlier experimental studies [48|, |49(, although attention 
has not focused upon these as pec ts therein. 

We have earlier suggested [501 ] that the vertical on-site potential anharmonicity, in combination with the high 
discreteness of dust crystals, may enhance energy localization via the formation of discrete breather excitations. The 
occurrence of DBs in ID dust crystals was investigated from first principles in Ref. [5l| . Our aim here is to extend those 
earlier results by elucidating the discrete dynamics of the transverse dust-lattice mode in 2D hexagonal crystals. A 
negative- valued "spring constant" is assumed where appropriate, in account of the inverse dispersion inherent in TDL 
vibrations. Our results are tested for values from real plasma discharge experiments. The analytical and numerical 
toolbox we used is described in full detail in earlier Refs. [I(J[l2>[H2] an d ^ s therefore only briefly summarized here. 

A number of comments are in row, to clarify our motivation and methodology. The scope of this article is two-fold. 
On one hand, we aim at presenting a generic investigation of the occurrence and stability of discrete lattice modes in 
hexagonal crystals, from first principles. On the other hand, the occurrence of discrete breather excitations associated 
with transverse dust lattice vibrations will be established through our work and lies in its very motivation. This 
double-sided scope is reflected in the structure of the text, which first adopts a generic hexagonal lattice formulation 
and then focuses closer on dust-lattices. For the sake of reference and rigor, regarding the former component (generic 
hexagonal lattice dynamics), we may cite a number of existing studies of hexagonal lattice configurations. Those earlier 
works have nevertheless focused on coupling (interaction) nonlinearity (e.g. via the Fermi-Pasta-Ulam paradigm, in 
fact disregarding on-site potential nonlinearity, a basic element in our case) [H, [54| or have employed the discrete 
nonlinear Schrodinger (DNLS) generic approach 1(| [l?], [H| • Of particular interest is the study of localized structures 
in hexagonal photonic lattices; e.g. [55l . l56l |57 |. The nonlinear Klein-Gordon approach that we adopt here was 
considered in Ref. [151 ] , upon which we have relied in our methodology. Nevertheless, we here generalize by considering 
the occurrence of breathers in systems wherein linear modes obey an inverse dispersion law; for this purpose, we have 
adopted a formulation which leaves the sign of the coupling coefficient arbitrary (i.e., cither positive or negative; 
read below). As regards the latter component we focus on (dusty plasma crystal transverse vibrations), earlier works 
[331 [36| have adopted a continuum approximation, whose validity may be questioned in certain real experimental 
situations involving highly discrete DP crystal configurations. 

The layout of the article proceeds as follows. Section [TT| is devoted to the analytical formulation of a general model 
for hexagonal lattice dynamics. The existence and stability of discrete breather type excitations in such a lattice is 
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investigated in Section HTT1 from a general point of view. Section lTVl is devoted to the analytical modelling of hexagonal 
dusty plasma crystals, in particular, while the generic results of Sec. IIIII are then explicitly applied to dusty plasma 
crystals in Section [V] Our results are finally summarized in the concluding Section. 



II. HAMILTONIAN DESCRIPTION OF A HEXAGONAL LATTICE 

Consider a hexagonal lattice like the one shown in Fig. [TJ Each site consists of a one-degree-of-frccdom Hamiltonian 
oscillator with coupled with its six nearest neighbors through a coupling constant e. 
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FIG. 1: The hexagonal lattice configuration is depicted, mounted with the site numbering employed in our model. The encircled 
oscillator is the "central" one. 



The Hamiltonian of the full system is 



i.j — — oo — — OO 



(4) 



where the i, j ' - indices denote the position of the oscillator in the lattice plane as it is shown in figHJ The variable 
Xij denotes the displacement of the particle located at site (i,j) in the direction vertical to the lattice plane and Pij is 
the corresponding momentum. A factor 1/2 is introduced in ([¥]) in order for each site pair entering the summation to 
be considered only once. Here it is assumed that the nonlinear on-site potential V of the oscillator possesses a stable 
equilibrium at (x,p) = (0, 0), with V"(0) = > 0. 



III. EXISTENCE AND STABILITY OF DISCRETE BREATHERS IN A HEXAGONAL LATTICE 



A. Existence and stability of single site breathers 



The existence of single site breathers can be investigated using the notion of the anticontinuous limit introduced in 
Q. This is constructed by taking e — > 0, i.e. by considering a chain of uncoupled oscillators. In this limit we consider 
one of the lattice oscillators, which we call "central" , moving in a periodic orbit of period Tt, while the rest of them 
lying at rest at (x,y) — (0,0). This state defines a trivially spatially localized and time periodic motion which is 
continued, for \e\ small enough, to provide a single site breather of period Tf,, if its period does not resonate with the 
phonon period, i.e. Ti, ^ mT p = 2n /uj p: \/m £ N (e.g. Q). The profile of a sample single site DB can be seen in figEl 

The linear stability of a breather is determined by the eigenvalues of the Floquet matrix (called also Floquet 
multipliers) of the corresponding periodic orbit, see e.g. 0, If all these multipliers lie on the unit circle of 
the complex plane then the periodic orbit is linearly stable, while if a multiplier lies outside the circle the orbit is 
unstable. Note that for every eigenvalue we also have its reciprocal and their complex conjugate(s), because of the 
Hamiltonian structure of the system (i.e. for every A^, there are also A" 1 , A* and A* -1 ). So, we cannot have just one 
eigenvalue outside the unit circle but only an even number of them. For e = 0, the above mentioned multipliers lie 
in two complex conjugate bundles at e ±JW P Tb , except of a pair of multipliers which lie at unity because of the phase 
degeneracy of the system. When the coupling switches on, for |e| ^ <C 1, the breather is formed. Correspondingly, 
the multipliers of the non-central oscillators move along the unit circle, without leaving it, because they are of the 
same Krein kind [59[, thus forming the phonon band, while the pair which lies at unity remains intact. We keep 



4 



increasing the value of |e| and the phonon band becomes wider. The DB, with the fixed period still exist and 
remains stable (changing thought its spatial profile) until the edge of the phonon band reaches +1. 

Our purpose is to investigate if, for the (experimentally) given potentials, the DB can exist for the values of e also 
provided by experiments. 



B. Existence of three-site breathers 

We shall now consider the anticontinuous limit, where for e — * we consider three "central" oscillators (marked by 
circles in Fig. [2]). These oscillators are moving in periodic orbits with the same period Tb = 2tt/u but in arbitrary 
phases, while the rest lie at rest. As it is proven [l5[, this state is continued for e 7 and |e| <C 1, if specific 
conditions for the phase difference between the central oscillators hold in the anticontinuous limit. The basic points 

i+3 M X X X X 

i+2 ~y V —*-y~~ --y- 

i+i — x — x — — K — x — 

\ / \ /J -1 \ 7 J \ /J +I \ / \ , 

i -V — ~-X- --X 
'' \ / j A / j \ / J +1 \ A +i \ / % 
i-i — X — X — X — x — x — 

\ / \ / j \ / \ , 

i-2 -X — X — X — X — X — X" 

FIG. 2: The hexagonal lattice, where the three "central" oscillators are shown. 



of the procedure which is fully described in [151 ] are shown below. 

Due to the phase degeneracy of the system, only the phase difference between the three central oscillators is 
significant in order to determine the specific periodic orbit for e = 0, which will be continued for e 7 to provide the 
three-site breather. This way, a natural canonical transformation to the central oscillators is induced 



u = wi 

h = w 2 

h = U> 3 



w 2 



A = J1 + J2 
h = J2 + J3 
h = J 3 



■h 



(•5) 



where we have defined the action-angle variables w, J. Note, that (pi denote the phase difference and since 03 = 
w 1 — W3 = —{4>i + (fe ), there are only two independent (j>iS, namely 4>\ and 02- Using these variables, and adopting 
the notation of Ref. (Hi , the condition for existence of three-site breathers becomes 



d(H x ) 



0. 



d 2 (H7 



70 



where 



(Hi) = f 



Hidt 



(6) 



(7) 



is the average value of H\ along the unperturbed periodic orbit. It is assumed that two conditions hold, namely: 
a) non-resonance of the breather frequency with the phonon frequency Xj, 7 mT p = m— and b) anharmonicity of 
single-site motion, viz. dw/dJ 7 0- 

The motion of a single uncoupled oscillator of the lattice in the neighborhood of its stable equilibrium point can be 
described by a cosine Fourier series in the form 



i(t) = ^ A w (Jj) COS 7 



TJ=0 



71=0 



A n (Ji) cos [n(uit + 



(8) 
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Since we consider the central oscillators to move with the same frequency u>b, the corresponding actions will have the 
same pertinent value Ji — J and so will have the Fourrier coefficients A n (Ji) = A n (J) = A n . Using ([5]), (Hi) becomes 

(Hi) = -- A n {cos(n0i) + cos(n0 2 ) + cos[n(0! + 4> 2 )}} , 



and condition becomes 



d(H ) 1 °° 

=oE nA « {sm(7i0i) + sin [n(0! + 2 )]} = 0. 



This is satisfied for all choices of harmonic content if Vn e N, 

sin(n0i) + sin \n{<j>\ + fo)] = 0, i = 1,2. 
At least three physically distinct solutions exist. The first one, with 

01,2 = 0, 



(9) 



corresponds to concerted motion of all three-sites oscillating in-phase, i.e. an in-phase three-site breather. The second, 
with 



01.2 = 7T, 

corresponds to an out of phase three-site breather. The final one, with 



n,2 



2tt 



(10) 



(11) 



describes periodic motion of all three-sites in a vortex like ring, to be henceforth referred to as a vortex breather. 
The spatial profile of these three configurations can be found in figs. 17191101 

C. Stability of three-site breathers 

As regards three-site breathers, for e = there are the two bundles of Floquet multipliers at e ±IWpTi> , which 
correspond to the non central oscillators, and in addition there are six Floquet multipliers at unity, a pair for each of 
the central oscillators. As it is already mentioned in the case of single site breathers, for |e| ^ the multipliers of the 
non-central oscillators spread along the unit circle, forming thus the phonon band. On the other hand, one pair of the 
multipliers that correspond to the central oscillators will always remain to unity, while the other two can move either 
along the unit circle or outside of it. If either one or both of these multiplier pairs leave the circle, the corresponding 
breather is unstable. On the other hand, if they both leave unity along the unit circle, the corresponding breather 
is linearly stable for values of |e| small enough. For increasing values of |e| the breather remains linearly stable until 
the multipliers of the central oscillators reach the phonon band. Then, since these multipliers are of different Krein 
kind, the breather is destabilized through a Hamiltonian Hopf bifurcation and a complex quadruple of multipliers is 
generated outside the unit circle. 

The Floquet multipliers of the central oscillators are given by = e aiT where oi are the characteristic exponents 
of the periodic orbit which correspond to the DB, and correspond to eigenvalues of the stability matrix E which in 
leading order of approximation is given by [l5j . 



E 
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i a 2 Hp , „ a 2 (gi) 

\ dhdh fc dl 2 dh 



b dtp 1 dl 2 
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a 2 (Hi) \ 

30i3<fe 

- £ -a^- 

c 9 2 (Hi) 

aii 302 

3 2 <Hi) 



(12) 



3/2 301 3/2 302 / 



If all the eigenvalues of matrix (|12|) lie on the imaginary axis, then the corresponding multipliers lie on the unit circle. 
If, in addition, these eigenvalues are simple to first order in e (or have definite "signature" , to be explained later) then 
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for small enough e the discrete breather is linearly stable. On the other hand, if any of the eigenvalue pairs of E is 
real then the corresponding DB is unstable. 

The various components of the stability matrix E are calculated as 



(13) 



with 



and 



1 oo 

/(</>) ^5> 2 ^cos(n0), (14) 



2 

n=l 



FHo _ to 
d 2 H I dJ 2 ~ dJ ' 3 



dlidlj 1 d 2 H du 



(15) 



J = * 



dJ 2 dJ ' 

The eigenvalues of E, i.e. the characteristic exponents of the DB, are to leading order of approximation, 



ffi,a = ±y [2/(20) + f(4>)] + 0(e), a 3A = ±yj -3e^ f(<f>) + O(e) . (16) 

We shall now investigate the behavior of these exponents for the 3 different modes introduced in the previous 
subsection. In doing so, we have to stress the significance of the product P = ejj (cf. the expressions above), which 
essentially determines the breather stability. We have chosen to keep the sign of the coupling strength e arbitrary, 
bearing in mind that inverse dispersive systems, such as transverse dust-lattice vibrations (to be discussed below), 
require for negative values of e to be considered in opposition to the usual Klein-Gordon configuration which is used 
for example for simulating the behavior of molecular crystals. 

a. Stability of in-phase breathers. The characteristic exponents for the first class of time-reversible solutions with 
0i = 4>2 = are to leading order of approximation 



*i,2,3,4 = ±y-3e^/(Q) + 0(e) , (17) 
Since 

/(0) = i5> 2 ^>0, 

71=1 

the sign of the product P — e^j of the coupling strength with the anharmonicity, determines the linear stability of 
the breather. If P < the specific solution is unstable, since all eigenvalues <Ji are real. On the other hand, for P > 0, 
the leading order calculation suggests linear stability (imaginary <7j). Because of the double multiplicity of <7j one 
could have expected that a splitting of the eigenvalues is possible, forming this way a complex quadruple which leads 
to instability. This, however, cannot happen due to symplectic signature reasons (the corresponding quadratic form 
is definite). This issue is thoroughly investigated in Refs. flBl [60|. 

Note that, the above stability results, as well as the ones that follow, apply both for e > (which corresponds to a 
classical nonlinear Klein-Gordon lattice) and for e < (which corresponds to Dusty Plasma Crystals). 

b. Stability of out-of-phase breathers. For the second class of solutions, corresponding to fa — n, the eigenvalues 
of the stability matrix E are 



-1.2 = - \j -erjj [2/(0) + /(tt)] + 0(e), a 3A = ±^-3^ /(tt) + 0{e). (18) 



Although the sign of the quantity 



2 

71=1 
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is not prescribed in general, one concludes from the above that 

2/(0) + /(tt)>0. (19) 

Therefore, in order to have both pairs of eigenvalues lying on the imaginary axis, i.e. for this class of solutions to be 
linearly stable, we shall require both conditions f(n) > and P = e^f > to be satisfied. 

Note that, in most cases we have f(ir) < 0, since the first term - which is usually dominant - is a negative term. 
Therefore, considering the exponential decay of Fourier coefficients of smooth functions, the rest of the terms usually 
do not change the negative sign of the sum. If this fact holds it means that the quantities f(n) and 2/(0) — f{n) will 
have opposite signs, so one of the Oi pairs will be real, which leads to instability. 

c. Stability of vortex breathers. The eigenvalues of the matrix E for the vortex configuration are 



<7i,2, 3 ,4 = ±^-3e^f/i +0(e), (20) 
with 

a-/(t)-'(t)-B- 1 )"- , ^-(t)- < 21 > 

Using the same arguments as in the previous case, usually f\ < 0. Then, if P > the corresponding vortex breather 
is unstable, while if P < the breather is linearly stable. 

IV. A MODEL FOR HEXAGONAL DUSTY PLASMA CRYSTALS 

Dusty plasma crystals form various configurations, the most common of which is a two-dimensional hexagonal 
lattice [2jj; see Fig. [TJ The typical inter-site spacing is of the order of (or slightly above) unity (in units of Debye 
length Xd)- Nearest neighbor interactions are therefore assumed in modelling dust crystals. 




FIG. 3: A snapshot of the DP crystal by a CCD camera. Courtesy of Dmitry Samsonov (University of Liverpool, UK). 

The vertical displacement z nm of the charged grain at site (i,j) in a DP crystal obeys an equation of motion in the 
form: 

-(F e -Mg) + F int , z -v-^ (22) 



We distinguish three contributions in the right-hand side (rhs), to be specified below. 

The first term corresponds to the the overall vertical force acting on a single grain (in the absence of coupling and 
damping) ; it consists of an upward force F e due to the electric field of the plasma sheath and of the force of gravity, 
Mg, pointing downwards. At the levitated equilibrium position of the crystal, at height zq, the two forces balance 
each other, viz. F e (zo) — Mg — 0. (We may set zq with no loss of generality, in the following.) Expanding the 
(variable) charge and electric field as q(z) w q + q' z + ^q'^(z)z 2 + ... and E(z) w Eq + E' z + \E'^{z)z 1 + ... (the 
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prime denotes differentiation with respect to z, while the index 
the overall vertical force reads 



! 0" denotes evaluation at the equilibrium position) 



F e (z)-Mg 



1 



Jq E' ( ; + 2q' () E' + q'iE () )z 2 

-Mui 2 (z + az 2 + bz 3 + ...) 
dV(z) 
dz ' 



(23) 



where the definition of all quantities is obvious. The transverse vibration eigenfrequency may be provided by a simple 
Bohm sheath theory [3(| [Hj]. The anharmonic character of the vertical on-site sheath potential V(z), which clearly 
appears in ab initio calculations for low density and low pressure [6l| and also in experiments (38l . [33 ] , is measured 
by the coefficients a and 6, which can be directly inferred experimentally [see (|25|) below; cf. Table 1]. The tilde was 
used in the latter expression to remind that these coefficients here have dimensions, to be henceforth removed via an 
appropriate scaling, viz. a — > roa>, /3 — > Tq/3. 

The vertical component Fi nt . z of the total interaction force Fi nt (due to the electrostatic coupling to neighboring 
sites) consists of the sum of the six nearest neighbor contributions, 



dU 
dr 

dU 
dr 



Szi 



,.=.tr'(' ,2+fc ?) ,/a 

r=a i=i i=i 



(24) 



where Szi = Zi — Zij (assumed -C a) and the dummy index I runs overs all six nearest neighbor - to (i,j) - indices in 
Figure f. For Debye interactions, viz. U(r) — e~ r / AD /r, the characteristic frequency ujt.q is obtained as in J5J). 

The last term in the rhs of (|22p represents energy dissipation (damping), and involves the collision frequency v. 
This term will be neglected, since v -C wo applies to number of real experiments [41| . The effect of damping will be 
included in future work. 

It is straighforward to combine the latter two expressions, (|23|) and (|2"4|) , into (|22|) . Furthermore, one may scale 

to obtain a reduced (dimensionless) 



space and time by the lattice spacing a and the inverse eigenfrequency ui Q , 



equation of motion. One thus obtains an equation of motion in the form M d Zij / dt 



the (classical) hamiltonian H is precisely given by setting Xi 



- dH (p i: j ,Zij)/ dzi j , where 



Pij 



MdZi 



^j/dt and lo p — 1 therein. Note 

however that the coupling coefficient in Q needs to be taken as negative- valued, viz. e = — w^q/wq. The negative 
sign of e, here arising naturally as an outcome of the algebra, was expected from the inverse dispersive character of 
vertical dust-lattice vibrations (see the negative slope in the observed dispersion curve, e.g., in Fig. 3 in [34|). 



DISCRETE BREATHERS IN A 2D DUSTY PLASMA CRYSTALS 



Consider a Hamiltonian of the form (H|) with e < and 

V{x) = x 2 /2 + ax 3 /3 + bx 4 /4: . 



(25) 



We shall consider three sets of typical values from real experiments [40l | for the nonlinearity parameters a, b and for 
the coupling constant e, shown in Table |TJ 





a 


b 


e 


set I 


0.01 


-0.04 


-0.034 


set II 


0.01 


-0.06 


-0.065 


set III 


-0.21 


-0.02 


-0.17 



TABLE I: Potential form coefficients occurring in laboratory experiments [40 



We have tested the numerical stability of the solutions obtained above, for the values of Table Q] Our numerical 
analysis is two-fold. It first focuses on the dynamics of the Floquet multipliers of the corresponding DB, for given 
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values of a, b and e, varying the latter in order to span a range of different values and thus identifying existing 
(in)stability regimes, in terms of the strength of the inter-grain electrostatic coupling. Once the stable/unstable 
regions were traced, a numerical simulation has been used to confirm our predictions for the breather stability. The 
main points of the analysis are presented below. 



A. Single site breathers 



We will first examine a system with the first set of values in Table fl] The on-site potential of the lattice has the 
form of Fig. |4]Ja) and the corresponding frequency range of an individual oscillator w(J) is shown in Fig. HJb). In 
order to compute a breather we have to consider one of the oscillators of the lattice to oscillate with a frequency 
ujb mojp. Note that in the specific case ui p — 1 due to normalization we have already performed. By taking a look in 
Fig. Hl^b) we can very easily choose a frequency Ub which does not resonate with the phonon frequency lu p . Although 
the graph Fig. S]Jb) is drawn with respect to action J there is a one-on-one correspondence with the amplitude of the 
oscillation which can be very easily computed even if we do not know the explicit form of the action-angle canonical 
transformation [521 ] . For example the results presented in Fig. [5] are for x max = 3.3 which results to tOb — 0.692. 




(a) 



(b) 



FIG. 4: (Color online) (a) The on-site potential of the lattice for set I of Table [I] and (b) the frequency ui of the oscillation for 
the same potential with respect to J. 



We can see that for |e| = 0.001 we have a pair of multipliers at unity while the rest lie in two very small conjugate 
bundles. For increasing values of |e| these bundles (which consist the phonon band) spread along the unit circle and 
finally for \e\ ~ 0.056 they reach +1. This means that the frequency of the breather has reached the phonon band, 
so the breather can no longer exist. So, we can expect the existence of single site breathers in such a system since 
the value of e in set I is \e\ = 0.034. Note that, in an infinite length lattice the phonon band would be continuous. In 
figs|5]it appears discrete because of the finite size of the system we use for computations. By using similar arguments 
we find that for sets II and III we can expect breather motions only up to |e| ~ 0.049 and e| ~ 0.037 respectively. 
This is somehow anticipated, since higher values of e the continuous behavior of the system. 




FIG. 5: 
in M. 



(Color online) A single site breather for |e| = 0.035. For a better representation, please refer to the video files provided 
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FIG. 6: (Color online) The distribution of the Floquet multipliers corresponding to the breather depicted in Fig. [5]for increasing 
values of lei. 



B. Three-site breathers 



We will consider again set I for our study. Due to EqfT7l and since < (Fig. 0Jb)) determines the linearly stable 
solution, at least for small values of e, is the in- phase configuration (Fig. [7|). 



(a) 



(b) 



FIG. 7: (Color online) An in-phase three-site breather for |e| = 0.015. For a better representation, please refer to the video 
files provided in [6^| . 

As it can be seen in Fig. [H the in-phase breather remains linearly stable until |e| ~ 0.017 when a Hopf bifurcation 
occurs and it destabilizes (Fig. [8]), which means that after approximately 120 periods of oscillation the energy starts 
flowing to the rest of the lattice. 




/ 0.5 




10 -0.5 


c a % 



|e| = 0.030 

FIG. 8: (Color online) The distribution of Floquet multipliers for increasing |e| fot the breather depicted in Fig. [7] The last 
two images show the destabilization of the breather for |e| > 0.017. 



Following the same kind of arguments we conclude that for the sets II and III of Table U the corresponding in-phase 
breathers remains stable until |e| ~ 0.016 |e| ~ 0.013 respectively. 

The out of phase (Fig. [9]) and vortex (Fig. [10]) configurations have [due to Eqs. (fT8j) . (|20| ] one and two pairs of 
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characteristic exponent which have real part to leading order of approximation, so they are linearly unstable for |e| 
arbitrary small. In Fig. [TTJ in can be seen that for |e| > the out of phase breather possesses a pair of multipliers in 
the real axis while the vortex breather possesses a complex quadruple. 




(a) (b) 

FIG. 9: (Color online) An out of phase three-site breather for e = 0.005. For a better representation, please refer to the video 
files provided in [6^| . 




(a) (b) (c) (d) 




(e) (f) (g) (h) 

FIG. 10: (Color online) A vortex breather for e = 0.008. For a better representation, please refer to the video files provided in 

VI. CONCLUSIONS 

In this article, we have established the occurrence in a hexagonal dusty plasma lattice of localized vibrational modes, 
in the form single-site or three-site breathers. Since transverse dust-lattice vibrations in a dusty plasma lattice are 
characterized by an inverse dispersion (backward wave), we have identified a phenomenological Klein-Gordon type 
Hamiltonian with negative coupling e and a quartic polynomial on-site potential to take into account this behavior. 
The parameters of the potential are given from real experiments and are divided into three sets originating from 
three different experimental configurations. We have considered the first set of values, since this best satisfies the 
discreteness requirement, having e — —0.034. For this configuration we conclude that single-site breathers can be 
supported, but the three-site ones are unstable for this value of the coupling. However, real values of the coupling 
parameter e in experiments span a large region of values, so that three-site breathers are expected to be stable for a 
different parameter set. 
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FIG. 11: (Color online) Floquet multipliers distribution for (a) the out of phase three-site breather depicted in Fig. [9] and (b) 
the vortex breather depicted in Fig. 1101 both for e = 0.008. Note that the both have multipliers outside the unit circle, even 
for small values of e, which imply instability. 
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